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Abstract. 

Recently the nonlocal coherent-potential approximation (NLCPA) has been 
introduced by Jarrell and Krishnamurthy for describing the electronic structure of 
substitutionally-disordercd systems. The NLCPA provides systematic corrections 
to the widely used coherent-potential approximation (CPA) whilst preserving 
the full symmetry of the underlying lattice. Here an analytical and systematic 
numerical study of the NLCPA is presented for a one-dimensional tight-binding 
model Hamiltonian, and comparisons with the embedded cluster method (ECM) 
and molecular coherent potential approximation (MCPA) are made. 
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1. Introduction 

For many years the coherent-potential approximation (CPA) pQ has been widely used 
for describing the electronic structure of substitutionally-disordered systems. However 
as a single-site mean-field theory, the CPA is not a fully satisfactory theory of disorder 
and it leaves much important physics out of consideration. In the CPA a site only 
feels the average effect of its environment and so, for example, statistical fluctuations 
in the chemical environment of a site responsible for sharp structure in the density of 
states (DOS) are not described 0. Furthermore the CPA is incapable of describing 
the effects of chemical short-range order (SRO) and other important environmental 
effects such as lattice displacements. However a new theory, the nonlocal coherent- 
potential approximation (NLCPA) [3], has recently been introduced by Jarrell and 
Krishnamurthy as a way of addressing such problems. The purpose of this paper is 
to investigate the general validity of the NLCPA by means of numerical calculations 
with a simple one-dimensional (ID) tight-binding model. 

There have been many previous attempts to generalize the CPA, though none 
have been completely satisfactory Of these theories, only the embedded cluster 
method (ECM) 0] and the molecular coherent-potential approximation (MCPA) [5] 
will be considered here. The ECM is a method for embedding a cluster of real impurity 
sites into a medium such at that provided by the CPA. Averaging over an ensemble of 
the various cluster configurations can therefore give a non self-consistent description of 
local environment effects. The MCPA goes a step further by enforcing self-consistency 
with respect to the cluster, and so yields a new effective medium. However this 
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medium has the periodicity of the cluster and so possesses the unsatisfactory property 
of breaking the symmetry of the underlying lattice. Such supercell periodicity also 
means the MCPA is computationally prohibitive since it involves a Brillouin zone 
(BZ) integration which scales with the cluster size, thus preventing its application to 
realistic systems. 

However the recently proposed NLCPA 3._ appears to possess all the attributes 
required for a satisfactory cluster generalization of the CPA. It was introduced in the 
context of a two-dimensional tight-binding model Hamiltonian as the static version 
of the dynamical cluster approximation (DCA) [7] used for describing short-range 
correlations in strongly-correlated electron systems. A recent application has been to 
describe impurity bound states in disordered d-wave superconductors . The NLCPA 
again determines an effective medium via the self-consistent embedding of a cluster, 
however the problem of maintaining translational invariance is solved by imposing 
Born- von Karman boundary conditions on the cluster, leading to an effective medium 
which has the site-to-site translational invariance of the underlying lattice. A further 
consequence is that the computational difficulties of the MCPA are alleviated since 
the BZ integration in the NLCPA does not scale as the cluster size increases. 

The NLCPA was subsequently derived by Rowlands [H] within the first-principles 
Korringa-Kohn-Rostoker (KKR) jl()l ITT] multiple scattering formalism, and the 
resulting KKR-NLCPA j5J method was first implemented for a realistic three- 
dimensional system (bcc Cu^Zn^) in IS. Moreover, the KKR-NLCPA method 
has recently been combined with density functional theory, enabling first-principles 
total energy calculations to be carried out as a function of SRO using self-consistently 
determined potentials Such self-consistent-field (SCF)-KKR-NLCPA calculations 
14 naturally include a description of the Madelung contribution to the total energy 
which is missing in the conventional single-site SCF-KKR-CPA, and important future 
developments include application to systems with spin |15j . strain |16| . and valency 
[T7] fluctuations. 

Due to the continued conceptual development of the NLCPA method and its first- 
principles application to realistic systems in particular, it is important to establish 
its general validity by addressing the following questions. Firstly, does the NLCPA 
produce physically meaningful results, in other words are the fluctuations described 
real? Secondly, does the DOS calculated using the NLCPA converge to the exact 
result in practice as the cluster size increases, and does it do so systematically? To 
date, these important questions have not been addressed. However this can be done 
by carrying out a systematic study of the NLCPA for a simple ID tight-binding model 
with diagonal disorder and nearest neighbour hopping. There are three reasons for 
this. The first is that in ID the exact result can be obtained against which meaningful 
comparisons can be made. Secondly, fluctuations are much more significant in ID 
and so detailed structure is expected in the DOS which can be accurately interpreted. 
Thirdly, such a model is computationally very simple and so large clusters can be 
considered, thus enabling the convergence properties of the theory to be investigated. 
Such an analysis is carried out in this paper, and any insight gained which could aid 
the interpretation of realistic calculations is presented. 

The outline of this paper is as follows. For completeness and for the purpose of 
introducing notation, the concept of the effective medium is described first in section[21 
Next, brief derivations of the CPA, ECM, and MCPA are given in sections EjTJ l2~2l 
and 12.31 respectively. In particular, the equivalence of the 'cavity' |S| and the 
'renormalized interactor' [2] formalism is established. An understanding of these 
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theories is crucial in order to appreciate and assess the advantages and possible 
limitations of the NLCPA. In section 12.41 a derivation of the NLCPA is given with 
an attempt to provide a physical picture of the theory, and comparisons with the 
formalism of sections 12.11 12.21 and 12.31 are made where necessary. In section |3 a 
numerical investigation of the NLCPA is presented, and results are compared with 
those of the CPA, ECM, and MCPA. Finally, conclusions are drawn in section 0] 

2. Formalism 

In a basis labelled by the sites of the direct lattice, consider a tight-binding 
Hamiltonian with matrix elements 

H ij =e i S ij + W ij (l-5 ij ) (1) 

and corresponding Dyson equation 

G lj = G% + G^t k G k] (2) 

k 

where e l is the on-site energy at site i, and W l ° is the parameter describing the 
hopping between site i and site j. For the case of a random substitutionally-disordered 
binary alloy with only diagonal disorder, the set of {e z } vary from site to site in a 
random fashion, taking the value eA with probability c or es with probability 1 — c. 
Denoting a specific configuration of {e 1 } by 7, and the corresponding Green's function 
for that configuration by G 7 , we are interested in (G 7 ), the average over all possible 
configurations. To do that we iterate © in a Born series, average term by term, and 
resum in the form 

(G» ) = G# + J2 G o^ H < G 7 ) ( 3 ) 
k,l 

where E u is the exact self- energy. J Since © is also a Dyson equation, it may 
be interpreted as describing an effective medium in which the site-diagonal part of 
E ,J decribes an effective on-site energy, and the site off-diagonal part decribes an 
effective correction to the hopping. The quantities (G^) and E u both possess the full 
translational (and point-group) symmetry of the underlying lattice, and so (|3J) may 
be expressed in reciprocal space as 

(G(k)} = Go(k) + G (k)S(k){G(k)) (4) 

The matrix elements of the Green's function are therefore given by the BZ integral 

(<%) = 7T- / dk ((Go(k))- 1 -E(k))- 1 e ik ( R -^) (5) 
"bz Jn BZ 

with (Go(k)) -1 = E — VF(k). Since it is not feasible to solve the problem exactly, the 
idea of an effective medium theory is to determine the best possible approximation to 
S« and (&/). 

I Clearly S y is not diagonal in the site index otherwise we would be making a mean-field 
approximation by writing the average of a product (e'G^J in terms of the product of the averages 



Investigation of the nonlocal coherent-potential approximation 



4 



2.1. Coherent-potential approximation 

The main approximation made by the CPA 1 is to assume a site-diagonal 
translationally-invariant self-energy 0e l Sij, where 0e* = 0e for all i. The CPA effective 
medium is then described by the equation 

0G y ' = G% + ^ G l o k 0e k 0G k: > (6) 

k 

In order to determine the medium, let us consider any site i. By removing the sum 
over all sites k and making up for the neglected terms by replacing the free particle 
Green's function with the cavity Green's function Q u , the site-diagonal part of © at 
site i can be formally rewritten in the form 

G" = Q H + g H 0e l 0G u = _1 - 0e 4 ) _1 (7) 

In fact Q" is related to a quantity known as the renormalized interactor A" [2] by 
0* = {E- A 11 )" 1 (8) 

where A" is given by the locator expansion |2] 

A" = W lk g k W kl + W lk g k W kl g l W u + ■■■ (9a) 

g* =(E- eT 1 (96) 

It can be seen that A" and therefore Q n depends only on the medium surrounding site i 
and is independent of the chemical occupation of i itself. It is therefore straightforward 
to define the Green's function for a real impurity embedded in the medium simply by 
replacing the effective site energy 0e l with a real site energy e l a at site i, where a = A 
or B. From (JZJ) this is given by 

G« = ((a")" 1 - el)' 1 = ( 0£ * - ej, + ^G'TT 1 ( 10 ) 
Now the CPA demands that 

;p a G£ = 0G" (ii) 



where P a is the probability that site i is of chemical type a. In other words, the 
replacement of an effective site energy by a real site energy should, on the average, 
produce no change to the CPA medium. Since the medium is translationally-invariant, 
it follows from JHJ that 0G" must also be given by the BZ integral 

G" = — !— / dk(E-0e-W(k)y 1 (12) 
Mbz Jn BZ 

The CPA medium is therefore determined from a self-consistent solution of (|11I12|1 . 



2.2. Embedded cluster method 

The ECM 4 is a method for embedding a cluster of real impurity sites in a medium 
such as that provided by the CPA. First consider © describing the CPA medium, but 
with the sites i,j lying within some cluster G. By restricting the sum over all sites k 
to run over the sites belonging to the cluster only, and making up for the neglected 
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terms by replacing the free particle Green's function with the cavity Green's function 
Q, @) can be formally rewritten in the form 

0G IJ = g IJ + J2g ik ^ k 0G kj (13) 

K 

where notation has been introduced so that all sites belonging to the cluster have been 
denoted by capital letters. In matrix form can be expressed as 

0Q CC = gCC + gCC 0s C 0a CC = ((gCCyl _ 0£ c)-i (14) 

where all matrices are in the space of the sites belonging to the cluster. In analogy to 
||SJ in the derivation of the CPA, Q_ cc is related to the cluster renormalized interactor 
A cc [2] by the expression 

g CG = (E- A CG )~ X (15) 

where A IJ , the matrix elements of A cc , are given by the expansion 

A IJ = Y W Ik g k W kJ + Y W Ik g k W kl g l W u + ■■■ (16) 

k<£C kgC,l(£C 

and g i is given by (|9&jl . It can be seen that A and therefore Q_ cc depends only on 
the medium surrounding the cluster C and is independent of the chemical occupation 
of C itself. The Green's function for an impurity cluster embedded in the CPA medium 
can now be straightforwardly defined simply by replacing the effective CPA cluster-site 
energy matrix 0e c with a matrix e c containing a configuration 7 of real site energies 
{e^} at each cluster site /, where a — A or B. From (|14(l this is given by 

G° C = ((G^y'-^y 1 = ( e c ~e c + (0G cc )- i y 1 (17) 

Once the CPA medium is known, the Green's function for an impurity cluster 
embedded in the medium can be straightforwardly determined through Q17[). and 
properties can then be calculated at any desired cluster site. However an average over 
all possible cluster configurations 7 is usually taken, in which case it will be found 
that {G^) ^ 0G^ for a cluster size iV c > 1 . This is due to an average effect on the 
electron hopping within the cluster caused by the disorder configurations. Thus 0e 
gains an off-diagonal correction and could be interpreted as a (non-self-consistently 
determined) cluster self-energy E c with matrix elements E IJ However this is defined 
only at the cluster sites, with the single-site CPA self-energy 0e l situated at all other 
sites, and hence there is no new effective medium. To do that, self-consistency with 
respect to the cluster must be enforced. One way forward is the MCPA jS], as described 
in the next subsection. 

2.3. Molecular coherent-potential approximation 

As described above, a cluster self-energy may be defined by averaging over an ensemble 
of impurity cluster configurations embedded in a medium. The idea of the MCPA [S] 
is to determine an effective medium with a periodically-repeating cluster self-energy. 
Such a medium may be defined by first dividing the lattice up into a collection of 
identical non-overlapping clusters, and then writing the Green's function in the form 

G CC ' = Q2 C ' + Y, GS°" (Z C " + W_ c " c ") 0G C " C ' (18) 

C" 
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These cluster matrices have site matrix elements 



0G CC ' 



G% c ' =G%, ieC,j€C' (19) 
j ij 

[£ C fc^].. = S«, [W^Scc]^ =W ij , i,jeC (20) 

The main approximation made here has been to assume that the (as of yet 
undetermined) cluster self-energy is cluster-diagonal. This means there are self-energy 
terms relating sites in the same cluster, but none relating sites in different clusters. 
The single-site translational invariance of the underlying lattice is therefore broken. 
However, the MCPA does require the medium to be invariant upon translation by a 
cluster so that it has the periodicity of a 'super-lattice'. As such, E c is independent 
of the cluster label G, and the Green's function 0G depends only on the difference 
between G and C. The real-space cluster dependence of Ijl8|l may therefore be 
removed using the Fourier transform 

0G CC (q) = 0G CC 'e-^ Rc - R c ) (21) 
c 

where Rc — R-c is the vector distance between the centres of clusters G and C', and 
q is a vector in the (smaller) BZ of the super-lattice. Equation (|18|l now becomes 

0G CC (q) = G CC (q) + G CC (q) (E C + W CC ) 0G CC (q) (22) 

where all matrices have the dimension of the cluster size only, and the inter-cluster 
dependence is expressed through the reciprocal-space vector q. 

The next step is to embed an ensemble of impurity cluster configurations into 
the medium defined above. This may done using the ECM of section l2~2l yielding the 
impurity cluster Green's function 

G^ c ^(^ c ~^ + ( G cc )- 1 y 1 (23) 

However, there are two main differences. Firstly, the cluster self-energy YF replaces 
the effective CPA cluster-site energy matrix e c . Secondly, the medium outside the 
cluster comprises of a periodically repeating cluster self-energy rather than single-site 
CPA effective site energies, and hence A cc has a slightly different expansion [2]. The 
final step is to generalize the CPA argument and demand that the average of G^ c 
over all 7 be equal to the Green's function for the medium itself i.e. 

{G? t C ) = 0G CC or P i G " = 0G " V/, J G G (24) 

7 

For a cluster containing N c sites, the number of configurations will be 2 for a binary 
alloy. Since the medium has the periodicity of the cluster, it follows from 12 11221) that 
0G must also satisfy 

G cc = J-[ dq(£-£ c -W cc -W CC (q)) _1 (25) 
where the integral is over the first BZ of the super-lattice. Therefore the MCPA 



medium is determined from a self-consistent solution of (I24l2f: 

The MCPA appears to be a natural generalization of the single-site CPA, but it 
does have some major disadvantages. Specifically, the intra-cluster hopping is now 
different from the inter-cluster hopping expressed through W_ cc (q) and so measuring 
properties at sites situated at the boundary of the cluster will in general give different 
results than sites at the centre of the cluster, for example. The supercell periodicity 
also causes problems in obtaining a well-defined spectral function, and leads to the 
computationally-demanding BZ integration in (|25|l which scales as N c increases. 
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Figure 1. (a) Real space ID tile (denoted by double-headed arrow of length 
4a) for a N c = 4 cluster. Sites are denoted by open circles, and a is the lattice 
constant, (b) Set of cluster momenta (periodic) denoted by closed circles for the 
N c = 4 cluster. The tiles centred at the cluster momenta are denoted by arrows 
and the solid line is the first BZ. The part of the tile centred at n/a that lies 
outside the first BZ can be translated by reciprocal lattice vectors into the first 
BZ to lie between —ir/a and —tv/Aa. (c) Set of cluster momenta (anti-periodic) 
denoted by closed circles for the N c = 4 cluster. Again the tiles centred at the 
cluster momenta are denoted by arrows and the solid line is the first BZ. 




2-4- Nonlocal coherent-potential approximation 

Like the MCPA, the NLCPA [2j determines an effective medium via the self-consistent 
embedding of a cluster. However, the key difference is that Born-von Karman 
boundary conditions are imposed on the cluster, leading to an effective medium which 
has the translational invariance of the underlying lattice. The first step is to solve the 
problem of an isolated cluster with Born-von Karman boundary conditions imposed. 

2.4-1- Cluster with Born-von Karman boundary conditions Conventionally in solid 
state theory the problem of a lattice comprising of a large number of sites with 
Born-von Karman boundary conditions is considered [18j . In ID, for example, this 
essentially means a very long chain of sites where one end of the chain maps round 
to the other end. If the lattice constant is a, then the first BZ extends from —ixja 
to +7r/a in reciprocal space, and the number of k points in the BZ is equal to the 
number of sites in the chain. It follows that if the length of the chain is now decreased 
to contain only say a cluster of N c sites, then the number of k points in the BZ will 
then equal N c . These are referred to as the set of cluster momenta {K n }, where 
n = 1, ..,N C . However, the boundaries of the BZ remain the same since the lattice 
constant is unchanged. The real-space cluster sites {/} (denoted by capital letters) 
and the corresponding set of cluster momenta {K„} satisfy the relation 

1 £ e iK„(R,-R, )=( ^ (26) 
C K„ 

for which the conventional lattice Fourier transform is recovered when N c — > oo. For 
example, in ID the (periodic) solutions are 

{R/} = (I-l)a 1=1,.., N c 
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Figure 2. (a) Self-energy E J (denoted by double-headed arrows) for the case of 
a N c = 2 cluster in the MCPA, for I ^ J. The vertical lines denote the boundaries 
of the tiles, (b) Same as (a) but for the NLCPA. 



for N c even. To date, only periodic solutions have been considered in the literature. 
However, for a given cluster size there is also a solution of corresponding to 
anti-periodic Born-von Karman boundary conditions (where in ID, for example, the 
wavefunction at one end of the chain maps back to minus the value at the other end) . 
In ID, the anti-periodic set of cluster momenta are given by 

. , P , (In - N c - 1\ 7T , , 

^=(nvr)a n=1 -^ (28) 

A ID example for N c — 4 is shown in figure 2] Notice that {K^} always includes 
the origin in reciprocal space, whilst {K^ p } are shifted (by -ir/(N c a) in ID) and 
lie symmetric about the origin. For the remainder of this derivation it is however 
convenient to simply choose one of the set of solutions above. A discussion concerning 
this choice is given in section !^. 4. 51 Cluster quantities that are translationally-invariant 
can be related in real and reciprocal space through i|26|) . for example for the cluster 
self-energy we have 



K„ 



E C ,(K„) = V E « e -«„<Hi-R.D (29) 



j 



2.4-2. Mapping cluster to lattice When mapping the cluster problem to the lattice, 
there are certain conditions on the geometry of the chosen cluster which must be 
satisified in order that that the translational and point-group symmetry of the lattice 
is preserved [3J- Specifically, if the cluster sites are now considered as 'embedded' 
in the lattice, then it must be possible to surround them with a tile which has this 
point-group symmetry and can be periodically repeated to fill out all space. This is 
analagous to the conventional Wigner-Seitz tile JS] used to surround a single site, and 
restricts the allowed values of N c for a given lattice. The construction for choosing 
appropriate clusters and tiles was first explained in [3] for a 2D square lattice and was 
generalized to realistic 3D lattices in § In ID, the tiles are simply lines (see 

figure 0) . 

Next, note that the periodicity of the real-space tiles define a reciprocal lattice, 
and each reciprocal lattice vector K will correspondingly be centred at a reciprocal- 
space tile. For N c = 1, this would correspond to the usual tiling of the conventional 
reciprocal lattice with Brillouin zones. However, for N c > 1, reciprocal space will be 
more densely filled with K points and the corresponding tiles will be smaller (by a 
factor of N c ). In the MCPA, such a tile would be the BZ of the 'superlattice' (see 

§ This construction would also apply to the MCPA if aiming to preserve point-group symmetry. 
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section . However, the idea of the NLCPA is to use N c such tiles to fill out the 
original BZ of the underlying lattice (see figure ^for a ID example), thus preserving 
translational invariance. In fact, the relevant K points will be the set of cluster 
momenta {K„}. 

The mapping of the cluster problem to the lattice may now be achieved as follows. 
In reciprocal space, consider the exact lattice self-energy E(k) of Q. The first step 
is to average E(k) over the momenta q within each of the N c tiles. This results in 
a 'coarse-grained' lattice self-energy E(K n ) which has a constant but different value 
within each tile. On the other hand, the cluster self-energy E C ;(K„) is defined only at 
the cluster momenta {K„}. The main approximation made by the NLCPA is to set 
E(K n ) to be equal to the value E c ;(K n ) within each tile n i.e. 

^— f dqE(K n + q) = E(K n )~E d (K„) (30) 

In other words, the lattice self-energy is approximated by that obtained from the 
cluster. The NLCPA approximation to the lattice Green's function in reciprocal space 
may now be represented by summing over the dispersion within each tile, yielding the 
set of coarse-grained values 

N I 

0G(K„) = -^/ dk(S-W(k)-E(K„))- x (31) 
"BZ Jn Kn 

which are straightforward to calculate since E(K„) is constant within each tile ^k„- 
Finally, using Ij26|l . the real space Green's function at the cluster sites becomes 

0G IJ = -^—Y[ dk(£7-W(k)-S(K n ))- 1 e ac - (Rl - R -' ) (32) 

The reason for coarse-graining the Green's function is to preserve causality by 
removing phase factors within the cluster, however these can be reintroduced after 
the medium has been determined (see section 12.4. 4fl . 

A physical picture of the resulting NLCPA medium may be seen by examining 
the consequences of approximating the exact self-energy E(k) by E(K„). Firstly, the 
self-energy is now restricted to act within the range || of the cluster sites {J, J} only 
(a more precise correlation length can be defined using Nyquist's sampling theorem 
|19l However, an electron can now propagate to any site in the lattice via the 
dispersion W(k), and will experience an identical self-energy E 7J at that site. Thus 
E remains a translationally-invariant quantity which depends only on the distance 
between sites / and J, now within the range of the cluster size, but independent of 
which site in the lattice is chosen to be site I. (E /J will also possess the point-group 
symmetry of the underlying lattice provided that the clusters and tiles have been 
chosen correctly). An illustration of this for the simplest case of a N c — 2 cluster in 
one dimension is shown in figure [3 

2.4-3. Impurity problem To determine the medium, the impurity problem must be 
solved. The first step is to define the reciprocal space cavity Green's function Q(K n ) 
via the Dyson equation 

0G(K„) = G(K n ) + £(K„)E(K n )0G(K n ) (33) 

| The self-energy only takes into account exactly nonlocal correlations that are within the range 
of the cluster size. However, when the medium is determined self-consistently, S IJ will include the 
effects of physics on a longer length scale at the mean-field level. 
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In diagrammatic terms Q(K. n ) is introduced to 'avoid over-counting self-energy 
diagrams on the cluster' [H]. Equation l|33f) can of course also be expressed in real 
space by applying the Fourier transform l|2tj[l to yield 

0G IJ = G IJ + J2 G IK Z KL 0G LJ (34) 

K,L 

where Q IJ is the real space cavity Green's function. In a similar manner to that 
introduced in the ECM and MCPA, here Q IJ is independent of the chemical occupation 
of the cluster itself. However, the important difference is that it is not possible to 
give an explicit real space expansion for Q IJ in terms of the cluster renormalized 
interactor as was possible in the ECM and MCPA. The reason for this is due to the 
translational invariance of the self-energy in the NLCPA, and is an example of the 
so-called 'embedding problem' 20 . If the ECM were used here to accomodate an 
impurity cluster in the NLCPA medium, then self-energy terms linking the cluster 
sites to the medium sites would have to be broken (see figure |5J). Such a situation 
does not occur in the MCPA since the self-energy is cluster diagonal and hence there 
are no self-energy terms linking sites in a cluster to the rest of the medium in that case. 
Although the NLCPA cavity Green's function Q IJ is a mathematically well-defined 
quantity, it is not the same quantity as that introduced in the ECM. The embedding of 
a cluster with Born-von Karman boundary conditions in the NLCPA medium should 
therefore be viewed as a mathematical construction rather than a true embedding in 
the conventional sense. Nevertheless, the NLCPA impurity cluster Green's function 
may be defined in analogy to H17|) by replacing the cluster self-energy in l|34|l with a 
particular configuration of site energies {e„} i.e. 

g!? = g IJ +Y. gIK ^ G i J ( 35 ) 

K 

The NLCPA self-consistency condition is equivalent to that of the MCPA given in 
|21 i.e. 

Y^P^lf = 0G IJ (36) 

7 

where P 7 is the probability of configuration 7 occuring. The effective medium 
is therefore determined from a self-consistent solution of (|32I36|1 . The cluster 
probabilities depend on the SRO parameter a and so SRO may be included by 
appropriately weighting the configurations in (|36() . provided that translational- 
invariance is preserved. Finally note that the NLCPA formalism reduces to the CPA 
for N c = 1 and becomes exact as iV c — > 00 since this would amount to solving the 
exact problem described by (|3I4(I . 

2.4-4- Calculation of Observables Once the medium has been determined through 
<|32I36[1 . there is no longer any need to coarse-grain the Green's function via l|31|l . 
Now, the Green's function may be calculated at any point in the BZ through 

0G(k) = (E- W(k) - S(K„)) _1 (37) 

and correspondingly at any sites i,j in the lattice by 

0G ij = _L I dk(E- W(k) - S(K„)) _1 e jk ( R '- R ^ (38) 
V-bz Jn BZ 
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In i|37l38[) above, E(K„) takes the appropriate value within each tile n. The 
configurationally-averaged DOS per site is given by the usual expression 

0n(E) = — Im 0G 11 (39) 

where 0G 11 may be calculated from (|32|) or from (|38|) above. However, when 
calculating site-off diagonal observables such as the spectral function, notice that 
S(K„) taking the appropriate (constant) value within each tile in i|37|) leads to 
discontinuities in 0G(k) at the tile boundaries. In order to remedy this, Maier et al 21 
proposed a scheme (within the DCA) to interpolate the self-energy using maximum 
entropy, however this can lead to causality being violated when the fluctuations are 
large. An alternative scheme for removing such discontinuities has been proposed by 
Batt et al |22 (see section 12 .4. 51 below') . 

2.4-5. Periodic and anti-periodic solutions It was mentioned in section [2.4.11 that 
for a given cluster size there are two possible sets of cluster momenta which may be 
used, {K p } and {K^ p }, corresponding to periodic or anti-periodic Born-von Karman 
boundary conditions respectively. In the literature to date, only periodic boundary 
conditions have been considered. The validity of calculations using only {K p } is 
analyzed in section^ However, it has been suggested that one should perform separate 
calculations using {K p } and {K^ p }, and then average over both results The 
philosophy behind this argument is that this has the effect of increasing the number 
of tiles, leading to a better sampling of the BZ. For example, in ID the number of tiles 
is effectively doubled, as can be seen from figure Q The validity of this suggestion is 
also examined in section |21 

However, Batt et al [22] have proposed a more sophisticated scheme for 
combining both solutions which also removes the discontinuities in 0G(k) mentioned 
in section 12.4.41 above. Firstly, let us label the Green's function (given by <|37|l ') for 
the periodic and anti-periodic solutions by 0Gp(k) and 0G^p(k) respectively. Next, 
observe that by construction 0Gp(k) is a better approximation to the exact result in 
the region of reciprocal space close to each of the points {K^} (in fact 0Gyip(k) has 
discontinuities at {K^}). Similarly 0GapO^) is a better approximation to the exact 
result close to each of the points {K^ p }. This suggests that, in brief, one should 
construct a new Green's function which follows 0Gp(k) close to each of the points 
{K p }, and follows 0G^p(k) close to each of the points {K^ p }, and interpolates in- 
between. Such a combined solution is guaranteed to be causal since it will always lie 
between the two causal extremes 0Gp(k) and 0G/ip(k). Furthermore, the problem of 
discontinuities is automatically removed since neither solution is followed where it has 
a discontinuity. Full details together with an implementation of this scheme are given 
in [22]. 

3. Results 

In this section results are presented for a ID tight-binding model with diagonal 
disorder only and nearest-neighbour hopping. Calculations for such a model were 
first performed with the NLCPA by Moradian |23] using N c = 4 and N c = 8 clusters. 
Here the random site energies can take the values €a = +2.0 and ep = —2.0, and 
the nearest neighbour hopping parameter is taken to be W = 1.0. Figure EH a ) and 
figure GUb) show the DOS results for a pure material comprising of energies ca or ep 
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Figure 3. (a) DOS (as a function of energy) for a pure material comprising of 
A sites, with tj\ = +2.0. (b) DOS for a pure material comprising of B sites, 
with es = —2.0. (c) Exact DOS results for a random A50B50 alloy of the pure 
materials above, (d) DOS for the same A50-B50 alloy obtained using the CPA. 



respectively. Notice that the pure bands just touch at E = 0.0, and so a theory for an 
alloy of the two pure materials must deal here with this difficult 'split-band' regime. 
Indeed, from the Saxon-Hutner theorem the gap at E = 0.0 for an alloy of the 
two pure materials should correspondingly be just vanishing. 

3.1. Exact Result and the CPA 

Figure 01c) shows exact DOS results for a random A 50 B 50 alloy of the two pure 
materials shown in figures IHIa),(b) obtained using the negative eigenvalue theorem 
|2fi| . The results have been displayed using centred histograms at 400 points from 
E = —4.0 to E = +4.0 i.e. an energy increment of 0.01 (a higher resolution yields 
too much structure for meaningful comparisons to be made). As expected, the DOS 
is highly structured and the gap at E = 0.0 is just vanishing. 

Next, figure E^d) shows a CPA calculation for the same A 50 B 50 alloy, displayed 
using centred histograms at the same energy resolution as the exact result. The 
CPA result is very smooth due to the neglect of nonlocal fluctuations in the disorder 
configurations. However, one of the central successes of the CPA is its ability to 
describe the split-band regime. Indeed the CPA correctly separates the DOS here 
into two sub-bands, but it is clear that it fails to correctly describe the onset of such 
split-band behaviour. This is because the contributions to the DOS near band edges 
are due to large clusters of like atoms which are not described by the CPA 0. 

3.2. Cluster theories 

The ability of the various cluster theories to improve upon the CPA result is now 
examined. All results here have been displayed using centred histograms at the same 
energy resolution as the exact result shown in figurc^c), and all energies have a 10~ 3 
imaginary part. 
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Figure 4. Configurationally-averaged DOS per site as a function of energy (in 
units of the bandwidth) for the various cluster theories (top to bottom) with 
cluster sizes N c = 2,4,6 (left to right). 




Figure 5. Configurationally-averaged DOS per site as a function of energy (in 
units of the bandwidth) for the various cluster theories (top to bottom) with 
cluster sizes N c = 8, 12, 16 (left to right). 
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Figures f4!5l show the DOS for the random A 5n B 5a alloy of figure using the 
various cluster theories (top to bottom) as a function of cluster size (left to right). 
For the ECM calculations shown in the first row, the configurational average has 
been taken over the 2 ° possible cluster configurations and the DOS measured on 
one of the central sites. The central site is likely to give a better description of the 
alloy since it is more influenced by the effects of the disorder than say a site on the 
boundary. Similarly, the MCPA calculations on the second row ('MCPA centre') have 
been measured at one of the central sites. For the MCPA calculations on the third row 
('MCPA average'), the DOS has been measured at each site on the cluster, and then 
an average taken over all such measurements. The fourth and fifth rows show NLCPA 
calculations using the periodic and anti-periodic sets of cluster momenta {K^} and 
{K^} respectively. Finally, the sixth row ('NLCPA average') shows an average of the 
periodic and anti-periodic results given in the fourth and fifth rows. 

First let us consider the ECM and MCPA central site results shown in the first two 
rows. As a general observation, it can be seen that all plots show increasing structure 
in the DOS as the cluster size increases. This arises from statistical fluctuations in 
the local environment of a site (smoothed out by the single-site mean-field CPA) due 
to specific impurity cluster configurations, and can be identified with corresponding 
structure in the exact DOS shown in figure E[c). Furthermore, the gap at E = 0.0 
is systematically filled in as N c increases due to larger cluster configurations present 
which contribute to the DOS at the band edges. Although this occurs more quickly 
with the MCPA than the ECM, it is interesting to observe that the ECM results and 
the MCPA central site results are very similar for all cluster sizes, indicating that 
self-consistency with respect to the cluster does not have a significantly large effect on 
the DOS results for this system (in the absence of chemical SRO). Indeed, the ECM 
and MCPA central site DOSs become identical for N c = 12, and both converge almost 
to the exact result at N c = 16. + Next, consider the 'MCPA average' results shown 
in the third row. For N c — 2, this is identical to the MCPA central site result by 
symmetry. As N c increases, more structure is again present in the DOS, although it 
is evident that the convergence towards the exact result is much slower. For example, 
the troughs in the DOS deepen slightly as N c increases from 8 to 16, but a higher 
cluster size would be needed for the troughs to reach the zero axis as seen in the exact 
result. 

Now let us turn our attention to the NLCPA calculations. First consider N c = 2. 
In the periodic (P) result, the peaks in the DOS can be identified with corresponding 
peaks present in the ECM and MCPA calculations, but the troughs present at E = ±3 
and E = ±1.5 are not seen elsewhere. The anti-periodic (AP) result is equivalent to the 
CPA here due to symmetry. On the other hand, the average over the P and AP results 
shown in the sixth row resembles the MCPA result much more closely, with peaks and 
troughs lying at the same energies. For N c — 4, both the P and AP calculations 
separately resemble the ECM and MCPA only over certain energy regions. However, 
the average over both the A and AP results again closely resembles those of the ECM 
and MCPA. Thus already, it is clear that P and AP calculations separately are not 
sufficient to yield a physical DOS, but the average does yield a physically meaningful 
DOS with peaks and troughs which can be associated with correlations in the cluster 

% An investigation of the ability of these cluster theories to deal with the effects of chemical short- 
range order is deferred to another publication. 

+ In theory, the MCPA has the incorrect boundary conditions as N c — » oo since the self-energy would 
become independent of k. In practice, however, the DOS converges quickly to the exact result. 
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disorder configurations. The reason behind this is clear from section 12.4.51 the P 
and AP calculations separately reproduce correlations arising from differing regions 
of reciprocal space, and thus a combination of the two calculations is necessary to 
reproduce all correlations within the range of the cluster size. The more sophisticated 
method for combining the A and AP solutions |221 mentioned in section 12. 4 .51 is likely 
to give even better results than the straight average displayed here, for example the 
'flat' regions in the N c = 2 calculation are likely to be smoothed out. Also, observe 
that as N c increases further, the A and AP plots appear increasingly similar, becoming 
almost identical at iV c = 12. This is expected since {K^} and {K^ p } become closer 
together as N c increases. 

Next, let us consider the convergence properties of the NLCPA. Importantly, 
it is clear that, provided the periodic and anti-periodic results are combined, the 
DOS calculated using the NLCPA converges smoothly and systematically towards the 
exact result. There are however clear differences in the convergence properties of the 
NLCPA compared to the ECM and MCPA. Firstly, the gap centred at E = 0.0 is 
filled in more quickly. Secondly, the convergence for the NLCPA is more systematic 
in that the maximum height of the DOS increases monotonically as N c increases 
(N c = 20 and N c = 24 calculations further confirm this). On the other hand, this 
means that large 'spikes' present in the ECM and MCPA DOS are not reproduced 
for small cluster sizes, although such spikes are seen in the exact result. Indeed, 
the ECM and MCPA plots show much more fine structure for a given cluster size, 
whilst the NLCPA results are much more 'smeared out'. Consequently, the ECM and 
MCPA central site calculations converge much more quickly to the exact result as N c 
increases. However, it is important to appreciate that the main (and more important) 
features of the exact result are reproduced fairly quickly by the NLCPA. Also, the 
ECM and MCPA calculations are very sensitive to the energy resolution at which they 
are displayed, unlike the NLCPA. These differences in the convergence properties are 
all a consequence of the fact that the NLCPA is based on Nyquist's sampling theorem 
[HUE]. Finally, observe that both the MCPA average (third row in figures EE} and the 
NLCPA have not yet converged to the exact result at N c = 16. Indeed there has been 
much debate in the literature |27U28ll2"§ll30| concerning whether the dynamical cluster 
approximation (DCA) jS] or the cellular dynamical mean field theory (CDMFT) 31 
(the analogues of the NLCPA and MCPA respectively for strongly-correlated electron 
systems) converges more quickly, carried out in terms of an average hybridization 
function coupling the cluster to the medium. Whilst a convergence comparison of the 
DOS using the NLCPA and MCPA average is not one of the objectives of this paper, 
it would be interesting to carry this out by using larger (and computationally very 
demanding) cluster sizes. It is clear, however, that a convergence comparison depends 
greatly on whether the quantity in question is sensitive to the energy resolution. For 
example, the convergence properties are likely to be much more similar for a quantity 
such as the integrated DOS. 

4. Conclusions 

A comparison of the formalism of the NLCPA with the ECM and MCPA has been 
presented. In the ECM and MCPA, the cavity function introduced has an explicit 
real-space expansion and so the coupling of an impurity cluster to the surrounding 
medium is well defined. It is found that this is not the case for the NLCPA cavity 
function (for N c > 1) as a consequence of the introduction of Born- von Karman 
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boundary conditions. This implies that the cavity should be viewed as a mathematical 
construction used to determine the medium, in departure from traditional effective 
medium theories. The consequences of introducing anti-periodic Born-von Karman 
boundary conditions within the NLCPA has also been investigated. Numerical DOS 
results for a ID model have been compared. It can be concluded that, firstly, 
the periodic and anti-periodic solutions must be combined when using the NLCPA, 
particularly for small cluster sizes. Importantly, only by then comparing the DOS 
with the exact result does it become evident that the NLCPA does indeed produce 
physical results and, for each cluster size, shares main features in common with both 
the ECM and MCPA. Since the NLCPA is based upon Nyquist's sampling theorem, 
it is however found that the DOS is far more 'smeared out' for a given cluster size. 
Whilst the NLCPA DOS does (importantly) converge systematically towards the exact 
result as the cluster size increases, it therefore does so less quickly than the ECM and 
MCPA (measured on the central site and in the absence of chemical SRO). On the 
other hand, such fine structure is less important and in any case would be smeared out 
in higher dimensions. More importantly, the NLCPA preserves the full symmetry of 
the underlying lattice, yields well-defined site-diagonal and site off-diagonal properties, 
and is computationally tractable for realistic systems. 
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